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ABSTRACT 

In  solidification,  when  the  process  is  limited  by  diffusion  of  released  latent  heat  or 
solute,  often  the  two-phase  interface  forms  finger-like  shapes  called  dendrites,  whose  tips 
are  nearly  paraboloidal  in  form.  For  a  pure  material  solidifying  into  an  undercooled  melt, 
if  surface  energy  and  gravity  are  negligible,  a  well-known  solution  due  to  Ivantsov  (1947) 
describes  the  steady  growth  with  a  paraboloidal  interface.  We  construct  a  regular  per¬ 
turbation  to  this  solution  for  a  downward  growing  axisymmetric  dendrite,  based  on  the 
smallness  of  a  buoyancy  parameter  G.  to  examine  the  effects  of  buoyant  flow  on  the  so¬ 
lidification.  The  analytic  solution  predicts  that  generally  the  buoyancy  enhances  growth 
and  changes  the  shape  of  the  interface,  giving  a  sharper  tip  and  wider  base.  These  effects 
depend  strongly  on  the  Prandtl  number,  and  also  on  the  Stefan  number  (undercooling). 
The  results  compare  well  v.ith  the  experiments  of  Huang  and  Giicksman  (1981a)  up  to 
G  ~  5000.  but  overpredict  convective  effects  for  higher  G. 


BUOYANT  CONVECTION  IN  DENDRITIC  SOLIDIFICATION 
INTRODUCTION 

Solidification  processes  are  important  in  a  variety  of  contexts,  and  are  of  particular 
interest  in  the  processing  of  materials.  The  small-scale  details  of  the  solidification  can 
affect  the  resulting  microstructure,  and  hence  the  bulk  properties,  of  the  final  product. 
This  practical  context  has  led  to  a  great  deal  of  research  toward  understanding  the  details 
of  solidification  processes. 

In  many  typical  solidification  situations,  an  initially  smooth  solid  surface  becomes 
morphologically  unstable  as  it  grows;  small  bumps  grow  into  long  fingers,  which  in  turn 
grow  side  branches,  yielding  a  tree-like  stucture  called  a  dendrite.  Experiments  have  shown 
that  the  growing  tip  of  the  dendrite  is  closely  approximated  by  an  axisymmetric  paraboloid 
for  many  materials  (or  an  elliptic  paraboloid  for  some  others). 

The  present  work  deals  with  the  effects  on  this  growing  dendrite  tip  of  buoyant  flow  in 
the  molten  phase,  driven  by  the  release  of  latent  heat  upon  solidification.  (With  the  pos¬ 
sibility  of  processing  materials  in  a  microgravity  environment,  prediction  of  gravitational 
effects  takes  on  practical  importance.)  To  gain  some  insight  into  the  effects  of  buoyancy 
on  dendritic  solidification  while  keeping  the  analysis  tractable,  the  simplest  situation  is 
considered;  we  analyze  an  isolated  axisymmetric  dendrite  growing  downward  into  an  un¬ 
dercooled,  pure  melt,  when  buoyantly  induced  velocities  are  small  compared  to  the  speed 
of  the  solidification  front  and  surface  energy  is  negligible.  We  find  a  solution  in  the  form 
of  a  regular  perturbation  expansion  in  the  small  gravitational  parameter  (G);  the  solution 
shows  explicitly  the  dependence  upon  the  Prandtl  number  ( P )  and  the  Stefan  number  (5, 
the  dimensionless  undercooling). 

There  has  been  much  previous  work  on  dendrites  with  convection,  although  a  rela¬ 
tively  small  proportion  applies  directly  to  natural  convection  near  isolated  dendrites.  The 
experimental  study  by  Huang  and  Glicksman  (1981a,b)  using  succinonitrile  show's  that 
buoyant  convection  can  significantly  affect  both  the  growth  speed  and  the  tip  radius  of 
isolated  dendrites,  when  the  growth  is  relatively  slow.  They  found  that  a  dendrite  growing 
downward  grew  faster,  with  a  smaller  tip,  than  one  growing  upward  in  the  same  condi¬ 
tions.  Similar  experiments  by  Tirmizi  and  Gill  (1987),  for  an  ice-water  system,  show's 
qualitatively  similar  results  with  platelike  dendrites,  and  with  dowm  and  up  switched  due 
to  water’s  increasing  density  with  temperature  near  freezing. 

Ananth  and  Gill  (1988)  present  a  theoretical  analysis  of  the  same  problem  considered 
here,  but  with  a  different  simplifying  approximation,  giving  a  different  region  of  appli¬ 
cability.  They  assume  the  dendrite  is  a  paraboloid,  and  essentially  apply  a  coordinate 
expansion,  formally  valid  within  a  fraction  of  a  tip  radius  from  the  very  tip.  They  de¬ 
rive  a  pair  of  coupled  nonlinear  ordinary  differential  equations  w'hich  include  the  nonlinear 
effects  of  buoyancy  in  this  region,  which  they  solve  numerically  for  some  choices  of  the 
parameters.  In  contrast,  we  assume  the  buoyancy  is  relatively  small  in  order  to  use  an 
expansion  in  the  gravitational  parameter  about  the  knowm  paraboloidal  solution  (Ivantsov, 
1947).  This  linearizes  the  buoyant  eonvection  effects,  which  allows  for  an  explicit  analytic 
solution  showing  all  parameter  dependence,  formally  valid  within  many  tip  radii  from  the 
very  tip.  This  solution  predicts  that  buoyancy  affects  not  only  the  growth  Peclet  number 
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(Pe,  the  dimensionless  product  of  growth  speed  and  dendrite  tip  size)  of  the  dendrite,  but 
also  its  shape;  the  perturbed  dendrite  is  not  paraboloidal. 

PROBLEM  STATEMENT 

Consider  a  smooth  dendrite  of  a  pure  substance,  growing  into  an  undercooled  quiescent 
melt.  We  assume  the  dendrite  is  axisymmetric,  with  the  axis  vertical,  growing  downward 
(see  Figure  1),  and  that  the  growth  is  steady  in  the  reference  frame  of  the  solid-liquid 
interface.  Surface  energy  and  non-equilibrium  effects  at  the  interface  are  neglected,  so  the 
liquid  solidifies  at  the  equilibrium  melting  temperature  Tm;  as  a  result,  the  entire  solid 
phase  is  isothermal,  at  Tm.  The  material  properties  in  each  phase  are  assumed  constant, 
with  no  density  change  upon  solidification.  The  Boussinesq  approximation  is  applied  in 
the  melt.  Then  the  governing  equations  in  the  melt  axe  the  conservation  of  energy,  mass, 
and  momentum: 


u  •  VT  =  kV2T  ,  (1) 

V  ■  u  =  0  ,  (2) 

(u  •  V)u  =  — -Vp  +  i/V2u  —  a(T  —  Too)g  ,  (3) 

P 

with  the  boundary  conditions 

at  the  interface:  T  —  Tm  ,  (4) 

u-Vz  ,  (5) 

—pL  Vz  ■  n  =  — fc(  VT)  •  n  ,  (6) 

deep  into  the  fluid:  T  —*  T oo  ,  (7) 

u  —y  Vz  ,  (8) 


where  T  is  the  temperature,  Tm  is  the  equilibrium  melting  temperature,  T ^  is  the  (under¬ 
cooled)  temperature  of  the  undisturbed  melt,  k  is  the  thermal  diffusivity,  k  is  the  thermal 
conductivity,  a  is  the  thermal  expansion  coefficient,  L  is  the  latent  heat  (per  unit  mass),  u 
is  the  velocity  (relative  to  the  interface),  V  is  the  axial  growth  speed  of  the  dendrite  (to  be 
determined  in  the  solution),  p  is  the  density,  p  is  the  reduced  pressure,  v  is  the  kinematic 
viscosity,  g  is  the  acceleration  of  gravity,  z  is  a  unit  vector  upward,  and  n  is  a  unit  vector 
normal  to  the  interface  (outward).  Condition  (6)  says  that  the  latent  heat  released  at  the 
interface  must  be  removed  by  conduction  into  the  melt,  while  (5)  and  (8)  say  that  the  solid 
and  the  undisturbed  melt  are  motionless  in  the  laboratory  reference  frame. 

To  nondimensionalize  the  problem,  a  length  or  velocity  scale  is  required.  Near  the 
dendrite  tip,  the  relevant  length  scales  are  the  tip  radius  P,  the  thermal  boundary  layer 
thickness  k/V,  and  the  viscous  boundary  layer  thickness  vjV.  All  of  these  involve  one 
of  the  unknowns,  either  tip  radius  R  or  growth  speed  V'.  As  a  result,  the  dimensionless 
solution  predicts  the  growth  Peclet  number  Pc  =  RV/ 2k,  but  cannot  predict  R  and  V 
separately.  This  problem,  of  selecting  which  of  a  family  of  solutions  is  realized  in  nature, 
is  inherent  in  the  problem  as  posed  (due  to  the  neglect  of  surface  energy),  and  is  common 


to  the  orginal  (non-buoyant)  solution  of  Ivantsov  (1947),  the  analysis  of  Ananth  and  Gill 
(1988),  and  the  present  work.  For  a  unique  physical  solution,  an  additional  selection 
criterion  is  needed,  e.g.,  based  on  the  morphological  stability  of  the  dendrite  tip. 

We  nondimensionalize  the  temperature  by  X  =  (T  —  T00)/(Tm  —  T^),  the  velocity  by 
u  =  u/F,  and  lengths  by  the  thermal  boundary  layer  thickness  x  =  x/(k/F),  then  drop 
the  tildes.  We  introduce  an  axisymmetric  stream  function  'P  to  conserve  mass,  and  use 
a  vorticity  formulation  to  eliminate  the  pressure.  Then  the  dimensionless  equations  for 
energy,  stream  function,  and  vorticity  are 


u  ■  VT  =  V2T  , 

(9) 

" = _v2  (7^)  • 

(10) 

p^rr 

(u-V)w-(w- V)u  =  PV2w-G— ^  , 

or 

(11) 

with  the  boundary  conditions 

at  the  interface:  T  =  1  , 

(12) 

$  =  r2/ 2  , 

(13) 

y\j>  —  r  r  , 

(14) 

z  n  =  5  (VT)  •  n  , 

(15) 

deep  into  the  fluid:  T  — »  0  , 

(16) 

>r2/ 2  , 

(17) 

where 

u  =  V  x  (  —<l>  \  ,  u  =  V  x  u  , 

(18) 

{r.  d,  z)  are  cylindrical  coordinates  with  corresponding  unit  vectors  r,  <f> ,  z,  and  the  three 
dimensionless  parameters  are  Stefan  number  (undercooling)  5  =  ( Tm  —T^Cp/L,  Prandtl 
number  P  =  vfn,  and  a  gravitational  parameter  G  =  ga(Tm  —  T^k/V 3  (like  a  Grashof 
number  times  Prandtl  number  squared),  and  Cp  is  the  specific  heat.  The  undercooling 
5  must  be  smaller  than  1  for  the  rate  of  solidification  to  be  limited  by  thermal  diffusion, 
as  assumed  here;  otherwise,  the  solidification  is  very  rapid  and  attachment  kinetics  limit 
growth. 

The  problem  formulated  above,  of  solidification  of  a  pure  substance  into  an  under- 
cooled  melt,  limited  by  diffusion  of  released  latent  heat,  is  analogous  to  a  solidification 
problem  for  a  binary  alloy,  under  certain  assumptions.  For  a  two-component  mixture,  the 
solid  and  liquid  phases  in  local  thermodynamic  equilibrium  at  the  interface  generally  have 
different  compositions.  For  example,  considering  the  major  component  the  solvent  and  the 
other  the  solute,  typically  the  solute  is  less  soluble  in  the  solid  phase,  i.e..  its  concentration 
is  smaller  in  the  solid  than  in  the  liquid.  Thus,  as  solidification  proceeds,  not  only  is  latent 
heat  released,  but  also  the  excess  solute.  When  the  mass  diffusivity  D  is  much  smaller 
than  the  thermal  diffusivity  k  (the  usual  case)  and  the  latent  heat  L  is  not  large,  then  the 
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speed  of  solidification  is  limited  by  diffusion  of  the  rejected  solute  into  the  melt,  and  the 
process  is  nearly  isothermal. 

Formally,  if  latent  heat  and  thermal  variations  are  negligible,  then  this  alloy  problem 
reduces  to  the  identical  dimensionless  form  above,  with  the  following  correspondences: 
T  —*  C,  k  —*  D,  S  (C,-C00)/(C,-C4),  P Sc  =  u/D,  G  -  g^Ct-C^D/V3,  where 
C  is  concentration  of  solute,  with  values  Cs  in  the  solid,  Cj  in  the  liquid  at  the  interface, 
and  Cqo  deep  into  the  melt,  D  is  the  diffusivity  of  the  solute  in  the  melt,  Sc  is  the  Schmidt 
number,  and  /?  is  a  solutal  coefficient  of  expansion,  defined  by  (p  —  Poo)/ Poo  =  0(C  —  C0 0). 
Therefore,  the  perturbation  solution  derived  below  applies  to  this  binary  solidification 
problem  controlled  by  mass  diffusion,  as  well  as  the  case  of  a  pine  substance. 

PERTURBATION  SOLUTION 

This  problem  presents  several  difficulties:  the  interface  between  the  solid  and  liquid 
phases  is  a  free  boundary,  the  energy  flux  condition  at  the  interface  (15)  is  nonlinear,  and 
the  convection  terms  in  the  transport  equations  (9)  and  (10)  are  nonlinear.  However,  when 
buoyancy  is  negligible  (G  — *  0),  a  solution  is  well  known,  due  to  Ivantsov  (1947).  In  this 
case,  there  is  no  motion  of  material  (in  the  reference  frame  of  the  quiescent  melt),  and 
the  solid-liquid  interface  is  a  paraboloid  of  revolution,  moving  steadily  along  its  axis  as 
the  liquid  solidifies.  In  the  reference  frame  of  the  interface,  the  flow  is  uniform  everywhere 
(u  =  z). 

To  get  some  indication  of  the  effects  of  gravity  while  avoiding  the  difficulties  of  the 
nonlinear  convection  terms,  we  consider  the  case  where  buoyancy  is  relatively  small  (G 
1),  and  seek  a  solution  as  a  small  perturbation  to  the  Ivantsov  solution. 

The  equations  suggest  a  regular  perturbation  in  G: 

T  =  T0  +  GT}+---  , 

^  =  ’J'o  +  +  •  •  •  ,  (19) 

u)  =  u>o  -f  G  u>i  +  ■  •  • 


Then  the  lowest  order  quantities  are  the  Ivantsov  solution,  which  is  most  easily  expressed 
in  paraboloidal  coordinates  (£,7,  <t>)  (see  Figure  2),  where 

Vr2  +  z2  +  z  Vr2  +  z2  —  z 

(  = - j -  ,  1  = - 5 -  .  (20) 

so  that  r  —  ,  2  =  £  —  77  , 

and  <t>  is  the  azimuthal  angle.  (For  comparison,  these  are  not  the  same  parabolic  coordinates 
as  in  Ananth  and  Gill,  1988.)  The  thermal  field  and  uniform  flow  are 

To  =  .  #0  =  j  =  2ft  ,  uo  =  0  ,  (21) 

and  the  position  of  the  interface  is  given  by 

rj  =  H0  ,  such  that:  H0  eWo-Ej(f/0)  =  S  ,  (22) 
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i.e.,  to  lowest  order  the  interface  is  a  paraboloid  whose  dimensionless  size  (or  Peclet  number: 
Pe  =  Ho  =  RV/2k,  R  is  the  tip  radius)  is  determined  by  the  undercooling  S.  (This  solution 
thus  represents  a  whole  family  of  physical  solutions  with  V  and  R  varying  inversely;  as 
mentioned  above,  an  additional  selection  criterion  is  needed  for  a  "nique  solution.) 

We  explicitly  consider  buoyant  effects  on  the  shape  of  the  dendrite  by  expanding  the 
interface  position  also  in  a  regular  perturbation  in  G: 

t1  =  H(O  =  H0  +  GH1(O  +  ---  ,  (23) 


Here  E\  is  a  standard  exponential  integral: 

e~x‘ 


ds 


At  the  next  order,  we  first  see  the  effects  of  buoyancy.  Linearizing  about  the  lowest 
order  solution,  the  equations,  in  parabolic  coordinates,  become: 


dTi  dT\  1  <9^i  , 


di 


dr]  2  di 


r,_l_  (  dTA  d(  m\ 

0  a?  di  )  +  a,  dn 


U\  = 


djj  j  duj  j 

t-ZT-V- 


~(Z  +  [  d^ 

P 


)  • 

'  a2*/ 

>  7?' 


dr\2 


di  '  drj 
with  the  boundary  conditions 

at  77  =  Hq  :  T\  +  H\  Tq  =  0  , 

T  5*!  „ 

*1  =  =  0  , 

dr] 

Hx+zm  =  -5 


HiU  +  H.HoT^  +  H0 


ar, 

J 


as  r]  — ►  oc 


T,  ->  0 


-»  0  , 


(24) 

(25) 

(26) 

(27) 

(28) 

(29) 

(30) 


where  primes  denote  ordinary  derivatives,  and  the  interface  conditions  have  been  trans¬ 
ferred  to  7/  =  Hq  by  expanding  T  and  in  Taylor  series  in  t]  about  H$. 

The  flow  at  this  order  is  driven  only  by  the  lowest  order  temperature  gradient.  So  we 
can  solve  for  the  perturbation  flow  first,  from  (25),  (26),  (28),  and  (30b),  then  use  the  flow 
to  find  the  temperature  and  interface  shape  perturbations. 

Because  the  driving  temperature  gradient  (Tq)  depends  only  on  one  coordinate  (77) 
and  the  interface  conditions  are  prescribed  on  a  coordinate  surface  (77  =  H0 ),  we  seek  a 
solution  that  is  separable  in  parabolic  coordinates.  Assuming  the  vorticity  is  of  the  form 


“AZ-rj) 


(31) 
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then  the  vorticity  equation  (26)  becomes 


The  general  solution  is 


Pf'+f 


-e-* 

riEi(Ho) 


(32) 


u>  1 


1 

(1  —  P)Ei(Hq) 


E2{t])  +  AE2{^)  +  B?) 


(33) 


where  A  and  B  are  arbitrary  constants,  and  E2  is  an  exponential  integral. 

The  corresponding  stream  function  is  not  in  simple  separable  form,  due  to  the  com¬ 
bination  (£  +  g)  in  the  denominator  of  the  coefficient  in  (25).  (That  factor  arises  from  the 
metric  coefficients  in  paraboloidal  coordinates.)  However,  assuming  a  stream  function  of 
the  form 

Vi(Z>rl)  =  Z29(v)  +  SHv)  ,  (34) 

then  (25)  becomes 

rig"  =  —2/  ,  .  (35) 

V  h"  =  -2(gf  +  g)  ,  (36) 


and  the  general  solution  is 

*> 1  =  f T^pkuh)  {eis'M  +  APgdP]  + 4 + ' Cu  +  D] 

HlhrW  4-  AP2hi(-~)  —  Cg2  +  2D^(1  -  log?)  +  Eg  +  F]}  .  (37) 

where  g\{g)  =  E2(g)  —  Ez(g),  h\{g)  =  2E$(g)  —  E2(g),  and  .4,  B ,  C,  D,  E,  and  F  are 
arbitrary  constants  (A  and  B  the  same  as  in  (33)). 

The  boundary  conditions  on  the  separable  parts  j(g)  and  h(g)  (i.e..  the  two  functions 
in  square  brackets  above)  are 


at  g  =  Ho  :  g  =  g  =  0  ,  h  =  h1  =  0  ,  (38) 

as  ?  — >  oo  :  g— >0  ,  h  — *  0  .  (39) 

Unfortunately,  no  choice  of  the  six  arbitrary  constants  can  satisfy  these  six  conditions, 
because  most  of  the  terms  do  not  decay  as  g  — >  oo.  Only  the  first  two  terms  in  each  set  of 
brackets  decay,  and  discarding  the  other  terms  to  satisfy  conditions  far  away  (39)  leaves 
only  one  free  constant,  A.  Conversely,  we  can  satisfy  no  slip  and  no  penetration  at  the 
interface  (38),  but  the  resulting  perturbation  flow  grows  without  bound  away  from  the 
interface.  Moreover,  whatever  solution  we  get  will  grow  without  bound  in  £,  i.e.,  along  the 
dendrite  far  from  the  tip.  (Looking  for  more  general  homogeneous  separable  solutions  to 
add  on  will  not  help;  they  will  not  have  the  right  asymptotic  form  as  £  — >  oc.) 
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Clearly,  a  solution  of  the  above  form  cannot  be  uniformly  valid  in  space.  However,  we 
can  find  a  locally  valid  solution  that  satisfies  the  interface  conditions.  This  solution  is  valid 
within  a  (dimensionless)  distance  of  0(G~l)  from  the  tip  in  all  directions;  we  call  this  the 
tip  solution.  (See  Figure  3.)  Along  the  dendrite,  feu-  from  the  tip,  the  nonlinear  effects  of 
buoyant  convection  become  important.  We  expect  a  strongly  convecting  boundary-layer 
in  this  region,  but  as  it  is  far  downstream  from  the  tip,  it  should  have  negligible  effects  in 
the  tip  region.  Far  from  the  dendrite,  a  distance  of  0(G-1),  is  well  outside  the  thermal 
and  viscous  boundary  layers  (unless  P  ~  G_1 ).  In  this  outer  region,  the  fluid  is  isothermal 
( T  =  0)  and  the  flow  is  irrotational;  it  matches  the  inner  flow  to  the  uniform  flow  of 
the  undisturbed  fluid.  That  is,  far  from  the  dendrite,  the  (potential)  flow  departs  from 
separable  form  to  conserve  mass.  So  we  seek  a  local  solution  about  the  tip  that  becomes 
potential  flow  far  away. 

We  choose  B  =  0  so  the  vorticitv  decays  in  77.  The  corresponding  terms  represent 
Pouiselle  pipe  flow,  without  the  pipe  (u  =  r2z).  We  also  choose  C  =  0  because  the 
corresponding  terms,  which  represent  stagnation  point  flow  (u  =  —  rr  -f  2zz),  grow  faster 
in  T)  than  the  remaining  terms.  Another  reason  to  discard  the  flow  terms  with  coefficients 
B  and  C  is  that  they  would  imply  that  the  dmdrite  induces  strong  flow  far  ahead  of  the 
tip,  which  is  physically  unreasonable. 

Then  the  other  constants  are  given  by  the  interface  conditions  (38): 

-g'i(Ho) 

‘  g[(H0/P)  ’ 

D  =  -[gl(H0)  +  APgl(H0/P)}  .  (40) 

E  =  ~[h\(H0 )  +  APh\(H0/P )  -  2D  log  H0]  , 

F  =  -{hl(H0)  +  AP2hl{H0/P)  +  2DH0{\  -\ogH0)  +  EH0}  =  0  . 

(It  can  be  shown  that  F  is  identically  zero.)  This  completes  the  tip  flow  solution. 

To  interpret  these  results,  note  that  the  flow  terms  <71(77 )  and  h\(rj)  in  (37)  repre¬ 
sent  buoyant  rotational  flow  driven  by  the  lowest  order  temperature  gradient:  they  decay 
exponentially  in  7;.  The  flow  terms  with  coefficient  A  describe  rotational  flow  (needed 
to  satisfy  110-slip  at  the  interface)  that  decays  exponentially  in  77/P,  i.e..  this  describes 
the  viscous  boundary  layer.  The  terms  with  coefficient  D  represent  a  mass  sink  dis¬ 
tributed  along  the  positive  z  axis  whose  strength  increases  linearly  with  r  (as  r  — »  0, 
u  —*  — (2r/r)r  +  log(r2/4z)z).  The  velocity  from  these  terms  grows  as  log 77,  i.e.,  far  away, 
the  dendrite  looks  like  a  mass  sink,  to  feed  the  rising,  growing  bound  urv  layer.  Finally,  the 
flow  proportional  to  E  is  uniform  flow  (u  =  z),  which  is  needed  to  satisfy  no  penetration 
at  the  interface.  So  outside  the  thermal  and  viscous  boundary  layers,  the  perturbation 
flow  becomes  irrotational,  as  expected. 

Using  this  flow  solution,  we  can  apply  the  same  approach  as  above  to  solve  the  thermal 
problem  (24),  (27),  (29),  and  (30a)  for  the  perturbations  to  the  temperature  (T))  and  to 
the  dendrite  shape  [H\ ).  Assume  solutions  of  the  form 

Ti(Ev)  =  tf-  i)p(i)  +  q(n) 

H,(Z)  =  HoU-DAo  +  A,} 


1 


(41) 


where  Aq  and  A\  are  constants  to  be  determined.  Then  the  problem  separates  to  give 


(Tip)' +  r)p  -p= -T^g  ,  (42) 

(<?«')'+  rn'  =  -Ti(g  +  ifc)  ,  (43) 

with  the  boundary  conditions 

at  Tj  =  Hq  :  p  =  -AqTq  ,  (44) 

q  =  -A,Ti  ,  (45) 

7A0  =  -S[-A0HQT^+p'}  ,  (46) 

A!=-S  C-Ao-A^HoT'o+q'  ,  (47) 

as  77  — >  oo  :  p  —y  0  ,  5  — >  0  ,  (48) 


where  9(77)  and  h(r])  are  the  flow  functions  defined  by  (34),  each  consisting  of  several  terms 
as  shown  in  (37),  and  the  identity  E[  +  rjE"  =  —riE[  was  used  to  simplify  the  boundary 
conditions. 

The  general  solution  can  be  found  by  applying  variation  of  parameters  for  each  of  the 
terms  in  g  and  h: 

-2 

P(v)  =  (i  P)E2(H0)  *Pl(- V)  +  APP^  +  dpM 

+  J(E2(ti)  -  EiM)  +  K[r,  +  1}}  ,  (49) 

q (7)  =  (2  _  p)e2(Ho)  +  Apq2(n)  +  DqM  +  Eq^iq)  +  lei(p)  +  m)  .  (50) 

where  J ,  I\,  L,  M  are  arbitrary  constants  and 
P\{q)  =  {4[£,2(2/7)  -  2Ei(2i7)]  -  ZEl{g)[E2(r])  -  £1(7)]  -  2e~'1[E2(q)  -  2E1{t])}}  . 

Pi(ri)  =  \{[E2(ri)  -  Ei(r!))[2PEi{~)  -  (1  +  2P)£3(^)  +  2 E2(^)  -  E^)} 

+  -p\p  +  1]  [G  22  ( ^7 )  —  G2\{t])  —  2Gl2(rj)  +  2Gu(t])  +  Go2(t 7)  —  Goi(7)]}  , 

P^Gl)  =  -{[£2(7)  -  £i(7)]U  -  log 7]  +  [7  +  1]£i(2)(7)}  , 

9i(7)  =  \E]( 7)  4-  e-'EM  -  2E\{2g)  ,  (51) 

72(7)  =  f?i(7)[d  -  P)E 3(^)  -  (2  -  P)Ej(^)  +  Ei(p 

+  p[(P  —  1)^21  (7)  +  (2  —  P)Gu{q)  —  Goi(t?)]  , 

<73(7)  =  2{£i(t/)[1  -  log 77]  +  e_r,(log 77  -  2)  -  E\2)( 7)}  , 

94(7)  =  -f-"  • 
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(52) 


and  E\  (77)  =  f  - ds  , 

J  v  s 

s 

Gmn(r})  =  J  Em(  —  )En(s)ds 

Applying  the  boundary  conditions  (44-48)  determines  the  constants.  For  the  temper¬ 
ature  to  decay  as  77  — >  00  then  K  =  M  =  0.  Manipulating  the  boundary  conditions  to 
eliminate  A0  and  A\  shows  that  J  and  L  are  determined  by 


HoP  (H0)  +  (2  +  Ho)p(Hq)  = 

H0q'{H0)  +  (l  +  H0)q(HQ)  = 

Then  the  shape  perturbation  is  given  by 

A0  =  Sp(H0)/H0 
A\  =Sq(H0)/H0 

This  completes  the  solution  for  the  0(G)  perturbation. 

P  =  1  the  some  of  the  forms  above  are  undefined,  but  the  limit  P  — ►  1  is  well  defined.) 

The  complicated  form  of  the  solution  makes  its  interpretation  difficult.  Interpretation 
can  be  aided  by  considering  the  opposite  limits  of  small  or  large  Peclet  number  H 0,  where 
asymptotic  forms  can  be  used  for  all  the  transcendental  functions,  leading  to  simpler 
results  at  the  leading  order.  While  both  limits  give  useful  approximate  results,  the  limit 
of  small  Hq  is  a  singular  limit,  all  the  different  terms  becoming  unbounded  at  0,  and  the 
the  asymptotic  forms  are  rather  complicated. 

For  large  H 0,  applying  the  relations  given  in  the  Appendix  yields  the  following  asymp¬ 
totic  forms  for  the  main  results: 


0  , 

-p(Ho)  . 


(53) 


(54) 


(For  one  particular  Prandtl  number 


-do 

A\ 


4(P  +  1)P0 
1 


17P  -f  7 
’  8(P  +  1)2P02 
11P  +  9 


+  0(W-3) 

Ho 


4  (P  +  1)  8(P  +  1)2P0 


+  0(WT) 
Ho 


(55) 


In  this  limit,  at  leading  order  .40  decays  in  Ho  and  Ai  approaches  a  constant  (maximum 
0.25). 

For  small  H0.  the  asymptotic  forms  are  more  complicated  due  to  logarithmic  singular¬ 
ities  at  the  origin.  The  results  are  given  as  quotients,  because  inverting  the  denominators 
would  give  an  infinite  number  of  terms  like  0(log-n(Po ))  before  terms  of  0(H 0). 
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A„  =  2  jlog(fl0)  |  log(2)  -  ^dilog(P  +  1)  +  log(P  +  1)  J  -  \  -  j|p) 

P,  ,^f~P2  ,  P'\  P *!  P  1] 

—log^J+logfP)^—  +  Ij__  +  ___  +  _j 
+  -jlog(2)log(P)  +  ?log(2)(7  +  1)  -  jdilog(P  +  1)(7  +  1) 

+  iog(P  + 1)(7  + 1)  (x  -  I  -  l  -  I5p)  -  7>°«2(p>(7  +  1) 

+  log(P)((7  +  l)(-^-+^)-^-i)+(7  +  l)((-2P  +  l)^-^)  +  0(P„log2(/f„))}  j 

(P  -  1)  {— 21og3(tf0)  4  log2(i?0)[21og(P)  -  67  -  5]  +  log(Ho)  [log(P)(47  +  3)  -  672  -  107  -  3] 

+  [dog(P)  -  7  -  1)  (272  +  37)]  +  O(H0  log3(^0))}  (56a) 


-42  =  | log2 (H0)  -- log(2)  +  Pdilog(P  +  1)  +  log(P  +  1)  (■ 


-2  P2  P  5 

T“+  2  +2+6P 


p  /op2 

+  vlog2(P)  +  log(P)(7r 


+  (2f-l)g+|(P-D 


rg  0 

+  l°g(#o)  3  log(2)  log(P)  -  log(2)(87  4  13)  +  Pdilog(P  4  1) (27  4  4) 

+  log(P  4  1)  4P7  4  474^-~4  2P+y4^)4  Plog2(P)  (7  4  2) 

+  Ios(P)  (j2  +  -2P~l)+  (2P  -  !)(T  +  2)y  4  (P  -  1)  4 


+  log(2)log(P)  (  —  46J  4log(2) 


-  2/2  -  j)  +  VP  -  !)(•>  +  2>y  +  (P  ~  1)  (y  +  f 

'-872  267  \ 

—3 - 3 - 6j  +  Pdilog(P  4  1)  (72  +  47  4  3) 


4  log(P+  1) 


±pli*  +  £4  +  272  +  52!  _  +  2p7  +  H2  +  I1-2£!  +  ?£  +  9  +  ± 

3  2  7  +  6  P  6  1  +  2  +  3P  2  +  2  +  2  +  2P 


+  ^log2(P)(72  +47  +  3) 

+  !og(P)  (  !l2  +  4  +  2£V  _  £4  +  H£4  _  2P7  -  2  +  2£!  -  44 

6'  ’  V  12  4  ^  3  2  6  1  6  2  2  / 

+(2P-l)(72  +47  +  3)^+{P-l)(^  +  H2  +  0]  +O(H0!og3(P0))|  / 

(P  -  1)  {21og4(P0)  4  log3(tf0)  [-21og(P)  +  87  4  7]  4  log2(P0)  [-  log(P)  (67  4  5)4  1272  +  217  +  8] 
4  log(tfo)  [-  log(P)  (672  4  IO7  4  3)  4  S73  4  217z  +  I67  4  3] 

4  [(-  log(P)  4  74  1)  (27  3  4  572  4  37)]  4  0  (H0  log  4(P0))}  (56b) 


10 


where  7  is  Euler’s  constant  (0.5772...),  and  the  dilogarithm  function  is  defined  by 

dilog(x)  =  f  l^±ds  . 

Jx  s  -  1 

following  Abramowitz  and  Stegun  (1972).  (The  lengthy  algebraic  calculations  preceding 
these  results  were  verified  using  a  symbolic  mathematics  program  called  Mathematica.) 

RESULTS 

The  above  perturbation  solution  describes  in  detail  the  buoyant  flow  and  the  resulting 
thermal  field  throughout  the  melt.  An  example  of  the  buoyant  flow  is  illustrated  in  Figure 
4;  in  the  reference  frame  of  the  immobile  solid  material  of  the  dendrite,  this  is  an  instanta¬ 
neous  view  of  the  complete  (unsteady)  flow  as  the  interface  sweeps  through.  The  velocity 
vectors  show  the  formation  of  a  rising  buoyant  layer  near  the  dendrite,  which  satisfies  the 
no-slip  condition  at  the  solid  surface.  In  this  example  ( H0  =  1,  P  =  2)  the  temperature 
field  decays  over  a  dimensionless  length  scale  of  about  1;  the  dropping  buoyancy  gives  a 
local  maximum  in  the  vertical  velocity  at  about  that  distance.  Farther  away,  the  flow  is 
drawn  inward  to  feed  the  growing  buoyant  layer.  Directly  ahead  of  the  dendrite  tip  there 
is  a  weak  upward  flow,  and  at  (horizontally)  farther  positions  the  flow  is  inward  and  even 
downward.  In  this  example,  the  net  effect  is  that  the  thermal  boundary  layer  (not  shown) 
is  slightly  compressed  at  the  tip,  leading  to  enhanced  tip  growth. 

The  main  results  of  interest  are  Ao  and  Ai,  which  describe  how  buoyancy  changes 
the  dendrite’s  growth  Peclet  number  Pe  and  shape.  The  dependence  of  A0  and  Aj  on 
the  parameters  is  shown  in  the  next  two  figures.  Figure  5  shows  the  dependence  on  the 
lowest  order  Peclet  number  H0 ,  which  is  is  monotonically  related  to  the  Stefan  number  5 
(dimensionless  undercooling)  by  the  Ivantsov  solution  (22).  The  perturbation  to  the  Peclet 
number  is  proportional  to  the  difference  of  the  two  curves: 

Pc  —  RV/2k  =  Jf|£=o  =  ffo[l  +  G{A\  —  Ao)]  ,  (57) 

where  G  is  the  buoyancy  parameter.  As  can  be  seen,  over  most  of  the  parameter  space 
( Hq.P ),  buoyancy  increases  the  Peclet  number  (Ai  >  Ao),  i.e.,  growth  is  enhanced  by 
convection  of  released  latent  heat  away  from  the  tip.  For  large  undercooling  (Ho  >~  1, 
S  >~  0.6)  the  difference  (Ai  —  Ao)  is  relatively  large,  though  in  practice  large  5  means 
rapid  solidification,  so  G  is  small. 

It  is  surprising  that  there  is  a  range  (P  =  0.2,  H0  ~  0.1)  where  the  curves  cross, 
predicting  a  very  slight  decrease  in  the  Peclet  number  from  buoyant  convection.  The 
interpretation  may  be  that,  when  viscosity  is  unimportant  (P  C  1)  and  the  dendrite  is 
moderately  small,  the  inward  and  downward  velocities  near  the  tip  compress  the  thermal 
boundary  layer  horizontally  more  than  the  weak  vertical  velocity  directly  ahead  of  the 
tip  compresses  it  vertically.  Then  the  net  effect  would  be  a  smaller  thermal  gradient  and 
slower  growth,  at  the  tip. 

The  departure  of  the  shape  of  the  dendrite  from  a  paraboloid  depends  only  on  Ao-  If 
it  were  a  paraboloid,  then  the  dimensionless  interface  position  f/(£)  would  be  constant,  so 
the  departure  is  H(()  —  Hq(  1  +  GA\ )  —  GHoAoH  —  1 ).  The  fact  that  Ao  is  positive  means 
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that  the  tip  (£  <  1)  is  relatively  smaller,  but  the  dendrite  widens  further  down  (£  >  1)  , 
because  the  buoyant  flow  redistributes  the  heat  away  from  the  tip  towards  the  base  of  the 
dendrite. 

Comparing  Figures  5a,  5b,  and  5c  shows  that  the  Prandtl  number  P  has  a  large  effect 
on  the  magnitudes  of  Aq  and  Ai,  as  well  as  some  effect  on  their  Ho -dependence.  Figure  6 
shows  the  dependence  on  Prandtl  number  P,  for  Ho  =  1,  i.e.,  S  =  0.596.  It  can  be  seen 
that  when  P  is  large,  buoyant  effects  are  relatively  small.  This  makes  sense,  in  that  the 
high  viscosity  inhibits  flow  in  the  thermal  layer,  near  the  no-slip  interface,  so  little  heat 
is  convected.  At  the  other  extreme,  the  results  become  nearly  independent  of  P,  because 
the  viscous  layer  near  the  interface  is  much  thinner  than  the  thermal  layer,  so  the  buoyant 
flow  is  limited  primarily  by  inertia. 

The  asymptotic  approximations  of  the  solution  for  large  and  small  H0  are  compared 
with  the  actual  solution  (for  P  =  2)  in  Figure  7.  (At  each  extreme,  both  the  leading-order 
approximation  and  the  next  higher-order  approximation  are  shown.)  It  is  apparent  that  the 
asymptotic  forms  converge  well,  so  are  useful  in  those  limits.  This  is  fortunate,  because 
the  actual  solution  becomes  difficult  to  calculate  at  both  extremes,  due  to  exponential 
dependence  for  large  H0  and  due  to  logarithmic  singularities  for  small  ff0. 

To  see  how  the  Ivantsov  solution  is  modified  by  buoyancy,  Figure  8  shows  growth 
speed  V  versus  tip  radius  P,  holding  constant  the  physical  parameters,  i.e.,  undercooling, 
material  properties,  and  gravity.  This  represents  the  family  of  possible  physical  solutions 
corresponding  to  one  particular  dimensionless  solution.  (The  lack  of  a  unique  physical 
solution  comes  from  neglecting  surface  energy  in  the  problem  statement.)  Then  the  form 
of  the  relation  of  R  and  V  is 

R  =  v  +  V*  ’  (58) 

where  A  =  2 kH0  and  B  =  2 ga(Tm  —  T<x>)k2 Hq(Ai  —  A0).  For  the  Ivantsov  solution, 
B  —  0,  giving  a  straight  line  on  the  logarithmic  plot,  and  for  large  V  the  buoyant  solution 
is  close  to  this  line.  But  for  slow  growth  (small  V,  large  buoyant  parameter  G)  the  buoyant 
solution  departs  from  the  Ivantsov  solution,  giving  relatively  faster  growth  for  a  given  tip 
radius  (or  a  smaller  tip  for  a  given  speed). 

Choosing  one  member  of  this  family  of  solutions  requires  some  selection  criterion 
outside  the  scope  of  this  paper.  One  criterion  that  seems  to  match  experimental  data  well 
(Huang  and  Glicksman,  1981a)  comes  from  considering  the  morphological  stability  of  the 
dendrite  tip;  the  criterion  is  of  the  form 


2kc/o 

V& 


=  a  (constant) 


(59) 


where  k  is  the  thermal  diffusivity,  d0  =  Tm~yCp/L 2  is  a  capillary  length,  Tm  is  the  melting 
temperature,  7  is  the  surface  energy,  Cp  the  volumetric  specific  heat,  and  L  the  volumetric 
latent  heat  of  fusion.  Different  theories  give  different  values  for  the  constant  <7,  but  the 
value  determined  experimentally  by  Huang  and  Glicksman  (1981a)  for  succinonitrile  is 
<7  =  0.0195.  This  selection  criterion  is  also  shown  on  Figure  8;  following  this  criterion 
shows  that  buoyant  convection  gives  both  faster  growth  and  a  smaller  tip  radius,  relative 
to  the  non- buoyant  (Ivantsov)  case. 


The  predictions  of  the  buoyant  theory,  augmented  by  the  selection  criterion  a  = 
0.0195,  are  compared  to  the  experiments  of  Huang  and  Glicksman  (1981a)  on  succinonitrile 
in  Figure  9.  The  values  used  for  the  physical  constants  are:  k  =  1.16  x  10-3cm2/s,  P  = 
22.8,  L  =  885.1cal/mole,  Cp  =  38.25cal/mole  K,  a  =  5.0  x  10-4K-1,  d0  =  2.7  x  10-7cm, 
and  g  =  980cm/s2  (from  Huang  and  Glicksman,  1981b  and  1981a). 

At  high  Stefan  number  5  (undercooling),  buoyant  convection  is  slow  compared  to  the 
interface  motion  (the  buoyancy  parameter  G  is  small),  and  so  buoyant  effects  are  negligible. 
As  the  Stefan  number  decreases,  the  growth  speed  V  also  decreases,  and  the  buoyant 
parameter  G  becomes  significant.  Both  theory  and  experiments  show  that,  relative  to  the 
non-buoyant  case,  buoyant  convection  increases  the  growth  Peclet  number  Pe,  enhances  the 
speed  V  of  growth,  and  decreases  the  radius  R  of  the  dendrite  tip.  However,  for  very  small 
undercooling,  the  theoretical  predictions  fail  to  match  the  experiments,  overpredicting 
convective  effects,  due  to  neglecting  the  nonlinear  inertial  terms.  Nonetheless,  considering 
that  the  buoyant  perturbation  solution  is  based  on  the  assumption  that  the  buoyancy 
parameter  is  small  (G  <  1),  the  range  of  agreement  is  surprisingly  good  (up  to  G  ~  5000). 

CONCLUSIONS 

Because  of  the  importance  of  solidification  in  the  processing  of  materials,  a  theoretical 
understanding  of  the  details  of  solidification  is  of  practical  interest.  Currently,  processing 
materials  in  a  microgravity  environment  is  possible,  but  very  expensive,  so  prediction 
of  gravitational  effects  in  solidification  may  be  helpful  in  deciding  when  the  expense  is 
justified. 

The  present  work  considers  one  special  case:  buoyant  convection  effects  on  an  iso¬ 
lated  dendrite  of  a  pure  substance,  solidifying  downward  into  an  undercooled  melt.  To 
further  simplify  the  problem,  it  is  assumed  that  the  dendrite  is  smooth  and  axisymmetric, 
that  surface  energy  effects  are  negligible,  and  that  buoyant  velocities  are  small  relative 
to  the  speed  of  the  solidification  interface.  The  dimensionless  problem  depends  on  three 
parameters:  Stefan  number  5  =  (Tm  —  T^Cp/L  (undercooling),  Prandtl  number  P,  and 
buoyancy  parameter  G  =  ga(Tm  —  T0 0)k/V3,  where  G  <C  1  by  the  last  assumption  above. 

The  buoyancy  parameter  G  is  the  product  of  a  Rayleigh  number  and  the  Prandtl 
number  P.  For  small  P,  inertia  limits  the  flow  and  G  gives  a  dimensionless  measure  of  the 
importance  of  buoyant  convection.  However,  for  large  Prandtl  number,  viscosity  limits  the 
flow,  and  a  better  measure  of  the  importance  of  buoyancy  is  the  Rayleigh  number,  G/ P. 
One  measure  that  covers  both  cases  is  G/(P  +  1). 

The  solution  is  constructed  as  a  regular  perturbation  in  G  about  the  non-buoyant 
paraboloidal  solution  due  to  Ivantsov  (1947).  The  first-order  (0(G2))  perturbation  quan¬ 
tities  describing  the  buoyant  flow  and  temperature  throughout  the  melt,  and  the  perturbed 
interface  position,  are  given  in  closed  form  above,  showing  explicitly  the  dependence  on 
S  and  P.  The  perturbation  solution  is  not  uniformly  valid  in  space;  it  applies  within  a 
dimensionless  distance  of  0(G-1)  1  from  the  dendrite  tip. 

The  main  result  of  interest  is  the  dimensionless  position  of  the  solid-liquid  interface 
in  axisymmetric  parabolic  coordinates  (£,77): 

interface  at:  V  =  Htf;  S,  P,  G)  =  H0(S)  {1  +  G[(£  -  1).40(5,  P)  +  A,(S.  P)]}  ,  (60) 
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where  Ho  is  the  Ivantsov  solution,  given  by  (22),  the  buoyant  perturbations  Ao  and  A\  are 
given  by  (54),  and  Figures  5  and  6  show  that  Ai  >  Ao  >  0  (with  the  exception  discussed 
in  the  Results  section). 

This  result  predicts  two  main  effects  of  buoyant  convection:  enhanced  growth  and  a 
modified  shape.  At  the  very  tip  of  the  dendrite,  £  =  0,  and  the  dimensionless  radius  of 
curvature  is  2 77,  so  the  growth  Peclet  number  is 

Pe=~=Hk=o  =  Ho{l  +  G(A1-A0)}  ,  (61) 

where  R  and  V  are  the  dimensional  tip  radius  and  growth  speed.  Thus,  because  G  and 
( Ai  —  Ao)  are  positive,  buoyancy  enhances  growth  for  a  downward-growing  dendrite.  (Fig¬ 
ure  5c  shows  the  exception:  for  small  P  and  a  narrow  range  of  Ho ,  the  difference  Aj  —  Ao  is 
slightly  negative,  predicting  slightly  reduced  growth  at  the  tip,  presumably  from  convection 
inward  and  downward  nearby.) 

The  dependence  of  the  dimensionless  interface  position  H  on  if  in  (60)  shows  that  the 
dendrite  is  not  a  paraboloid;  because  Ao  >  0,  the  tip  is  smaller,  and  toward  the  root  the 
dendrite  widens,  relative  to  a  paraboloid.  This  buoyant  modification  of  the  dendrite  shape 
is  a  unique  feature  of  the  current  work;  previous  analyses  assume  a  paraboloid  shape. 
Furthermore,  while  the  dimensionless  solution  cannot  predict  R  and  V  separately,  the 
use  of  any  selection  criterion  of  the  form  R2V  =  constant  (such  as  criteria  based  on  tip 
stability)  will  predict  both  faster  growth  and  a  smaller  tip  due  to  buoyant  convection  of 
released  latent  heat  away  from  the  tip,  toward  the  root. 

The  magnitude  of  these  buoyant  effects  depends  strongly  on  the  Prandtl  number  of  the 
melt,  as  shown  in  Figure  6.  For  large  P,  the  viscous  length  scale  is  larger  than  the  thermal 
layer,  and  the  no-slip  solid  interface  inhibits  convection,  so  buoyant  effects  are  relatively 
small.  However,  for  small  P,  the  viscous  layer  by  the  interface  is  thin  relative  to  the  thermal 
layer,  so  the  buoyant  convection  is  effectively  inviscid,  limited  by  inertia,  and  buoyant 
effects  axe  relatively  large.  Thus,  conclusions  based  on  experiments  or  calculations  for  a 
high-Prandtl  number  material  (e.g.,  succinonitrile:  P  =  23)  should  be  applied  cautiously 
to  the  processing  of  low-Prandtl  number  materials  (e.g.,  metals),  and  vice  versa. 

Due  to  the  neglect  of  nonlinear  inertial  terms,  the  present  work  overpredicts  buoyant 
effects  when  G  is  large.  However,  the  theory'  (with  the  selection  criterion  based  on  sta¬ 
bility)  matches  the  experiments  well  up  to  G  ~  5000,  despite  the  formal  assumption  that 
G  <  1.  In  comparison,  the  nonlinear  calculations  of  Ananth  and  Gill  (1988)  match  the 
experiments  well  even  for  large  G]  their  selection  criterion  involves  choosing  the  Grashof 
number  to  fit  the  experimental  Peclet  number.  Their  method,  which  assumes  the  dendrite 
is  a  paraboloid,  requires  numerically  solving  a  pair  of  coupled  nonlinear  ordinary  differ¬ 
ential  equations  for  each  set  of  parameters,  and  their  published  results  are  limited  to  one 
Prandtl  number  (P  =  23.1).  The  results  of  the  present  work  are  expressed  in  closed  form, 
and  are  relatively  easy  to  calculate,  involving  only  a  few  different  transcendental  functions 
and  several  arithmetic  combinations.  (Rough  predictions  can  be  made  by  interpolation 
using  Figures  5  and  6.) 

While  the  perturbation  solution  applies  only  for  a  downward  growing  dendrite,  the 
perturbations  (by  assumption)  are  linear  in  the  buoyancy  parameter  G.  and  it  is  tempting 
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to  apply  the  same  solution  to  an  upward  growing  dendrite  by  merely  reversing  the  sign 
of  G.  We  speculate  that  this  may  give  reasonable  results  if  the  magnitude  of  G  is  small 
enough,  so  that  the  speed  of  the  interface  as  it  sweeps  through,  freezing  the  material, 
overwhelms  the  instantaneous  local  fluid  velocity.  However,  for  larger  buoyancy,  i.e.,  slow 
solidification,  the  thermal  layer  of  fluid  rising  toward  the  tip  of  the  dendrite  could  keep 
rising  past  the  tip,  forming  a  plume  above  the  dendrite.  Such  a  plume  would  change  the 
character  of  the  flow  completely,  so  the  perturbation  solution  would  not  apply. 

The  buoyant  perturbation  solution  can  also  be  interpreted  to  apply  to  the  growth  of  a 
dendrite  of  a  binary  alloy,  if  the  process  is  controlled  by  the  mass  diffusion  of  solute  rejected 
at  the  solidification  front  (as  discussed  in  the  Problem  Statement).  This  would  be  the  case 
if  the  mass  diffusivity  D  of  the  solute  is  much  smaller  than  the  thermal  diffusivity  k  in 
the  melt,  and  if  the  latent  heat  is  negligible.  Then  the  temperature  is  effectively  constant 
(Too),  and  the  transport  problem  reduces  to  convection  and  diffusion  of  solute,  measured 
by  its  concentration  C .  To  interpret  the  solution  for  this  alloy  problem,  substitute  C  for 
T,  nondimensionalize  using  D  instead  of  k,  and  the  new  dimensionless  parameters  become 
S  ->  (Cj  -  Coo)/(C,  -  C3),  P  -4  Sc  =  u/D,  and  G  -  g0(C,  -  C^D/V3,  wh, C,  is 
the  concentration  in  the  solid,  C/  in  the  liquid  at  the  interface,  and  deep  into  the 
melt,  Sc  is  the  Schmidt  number,  and  0  is  a  solutal  coefficient  of  expansion,  defined  by 
(P  ~  Poo)/ Poo  =  P(C  -  Coo). 

The  methods  employed  here  may  also  apply  to  dendrites  shaped  more  like  elliptic 
paraboloids.  The  key  element  that  makes  the  perturbation  analysis  feasible  here  is  the 
existence  of  a  simple  solution  to  the  nonbuoyant  case,  depending  only  on  one  coordinate  in 
a  separable  coordinate  system,  that  gives  the  shape  of  the  free  boundary  and  accounts  for 
the  nonlinear  flux  condition  at  the  interface.  The  family  of  such  solutions  is  quite  extensive 
(Canright  and  Davis,  1989),  including  one  for  an  elliptic  paraboloid.  (In  fact,  solutions 
are  known  for  binary  alloys  when  both  thermal  and  solutal  variations  are  important,  or 
ternary  alloys  when  thermal  variations  are  negligible,  so  the  perturbation  approach  may 
be  applicable  to  such  cases  as  well.)  The  detailed  form  of  a  buoyant  perturbation  to  the 
elliptic-paraboloid  basic  state  would  be  more  complicated  than  the  solution  given  here,  in 
that  the  flow  would  depend  on  all  three  spatial  coordinates,  and  the  special  functions  may 
also  be  more  complicated.  We  speculate  that  the  qualitative  features  in  that  case  may  be 
similar  to  those  of  the  present  case. 
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APPENDIX:  SPECIAL  FUNCTIONS  IN  SOLUTION 

In  the  perturbation  solution,  the  various  functions  (g's,  h' s,  p’s,  q's)  describing  the 
stream  function  and  temperature  perturbations  are  expressed  in  terms  of  exponential  in¬ 
tegrals  En  and  integrals  thereof. 
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The  exponential  integrals  are  standard  functions,  and  can  be  found  in  standard  refer¬ 
ences,  e.g.,  Abramowitz  and  Stegun  (1972).  Some  of  their  properties  are  recounted  here, 
for  convenience.  The  definition  is 


;„<*)  =  J 


0°  c-« 


ds  , 


(62) 


where  we  are  only  concerned  with  non-negative  integer  values  for  n.  Some  special  cases 
are 


Eq(x)  = 


Ei( 


ds  . 


Differentiating  the  definition  shows  that  the  derivatives  are 

d 


dx 


En(x)  =  -En-i(x)  . 


Integrating  the  definition  by  parts  gives 

En+1(x)  =  -  [e-1  -  xEn(x)] 


(63) 


(64) 


(65) 


This  can  be  used  to  reduce  all  such  exponential  integrals  En  to  forms  involving  only  the 
special  function  E\.  At  the  origin,  E\  has  a  logarithmic  singularity 


OO  ,  ^ 

£i(i)  =  -iiogi+i-)  -  ’ 


(66) 


k- 1 

where  7  ~  0.5772  is  Euler’s  constant.  For  large  arguments,  the  asymptotic  series  is 

-*  —  k\ 


—  ,  ,  e  v — '  k: 

it =0  v  ’ 


(67) 


The  fmiction  E\  1  is  less  common,  but  has  applications  to  stellar  atmospheres.  It  was 
first  defined  by  van  de  Hulst  (1948);  we  use  his  notation: 


£'2,(x)  =  j‘ 


Ei  M 


ds 


Switching  the  order  of  integration,  this  can  also  be  expressed  as 


:(,2,(x)  =  J 


00  e—  log. 


■ds.  , 


(68) 


(69) 


which  is  more  convenient  for  numerical  calculation.  (Kourganoff,  1952,  Appendix  I)  shows 
this  result  and  many  others  related  to  exponential  integrals.)  Near  the  origin,  the  following 
expansion,  due  to  van  de  Hulst  (1948),  converges  quickly: 


E[2)(x)  =  ^(logx  +  7)2  +  ~  +  £ 


k=l 


i~*)k 

k2k\ 


(70) 
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For  large  arguments,  integration  by  parts  gives  the  following  asymptotic  series: 


E\2\x) 


(-I)*"1*!  A  1 

X*  j 

J=1 J 


(We  have  not  seen  this  or  subsequent  results  published  previously.) 

Van  de  Hulst  (1948)  defined  certain  functions  he  called  Gmn;  the  functions  we  call 
Gmn  are  similar  to,  but  different  from,  his  functions.  We  use  the  definition 


r°°  s 

Gmn(x)  =  J  Em(  —  )En(s)  ds  , 


where  the  parameter  P  is  the  Prandtl  number  for  our  purposes.  Following  the  methods 
used  by  Chandrasekhar  (1950,  Appendix  I),  i.e.,  integration  by  parts  and  the  identities  for 
En,  all  of  the  Gmn  we  need  can  be  reduced  to  simpler  functions: 

C02(x)  =  P2t~xlpEl{x )  +  (P  +  P2)Ei(ftx)  , 

G„(x)  =  e~z E\(x / P)  +  Pe~x/p Ei(x)  -  xE1(x)El(x/P)  -  (P  +  , 

G12(x)  =  \  \Pe~0x  +  (1  -  x)e~z E\(x / P)  -  (P2  +  Px)e~zl p E,{x) 

£ 

+  x2E1(x)E,(x/P)  +  {P2  -l)E,(Dx)}  , 

G2i(x)=j  e-^  -  ^{1+  x)e-E1(x/P)  +(P  -  x)e-/pE,{x)  (73) 

+  yE1(x)El(x/P)-(P-^)El(0x)  , 

Gjxfx)  =  1  (P  +  1  -  -  i(l  +  x  -  x2)e~rEi(x/P)— 

-  ( P 2  +  Px-  x2)e~'",E,(x)  -  yE1(x)E,(x/P)  +  ( P 2  +  l)£,(3x)  , 

where  ft  =  1  +  -p.  The  one  exception  is  G 01,  which  apparently  cannot  be  reduced  to 
standard  functions.  By  definition, 


—  Gq\{x)  -  J 


OC  -,/P 


■Ei(s)ds  . 


Switching  the  order  of  integration  yields  a  form  more  amenable  to  numerical  computation: 


JGnlx)~L  Vlog  (;-?)* 
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Near  the  origin,  the  series  representation  is 


■^Goi(x)  =  ^(logx  +7)2  -  -  log2  P  ~  dilog(P  +  1) 


T  1  T* 

i  //I _ i  i  /  d2  i  a  D  1 


+  (-(log!  +  7)  -  P  +  1)?  +  ((log!  +  7)  +  j(P2  +  iP  -  l))jpj  l76) 

+  (-(log  1  +  7)  -  g(2P3  +  OP2  +  18P  -  1))  jf^J  +  0(x4  log  i)  , 

where  the  dilogarithm  function  is  defined  by  dilog(t)  =  —  f*  following  Abramowitz 

and  Stegun  (1972).  (The  general  form  for  the  nth  term  is  not  apparent.)  Integration  by 
parts  gives  the  asymptotic  series  for  large  arguments: 


1  r  (  \  Pe  /3r  V' 
•pGcfx).  —  Z 

*=1 


y _ l _ 


Using  the  reductions  above,  all  the  perturbation  functions  (g' s,  h' s,  p’s,  q's)  can  be 
expressed  in  forms  where  the  only  special  functions  are  E\,  E\  ,  and  Goi: 

0i(x)  =  \  [(1  +  x)e~z  -(2x  +  x2).Ei(x)]  , 

hi(x)  =  —  xe~z  +  (x  +  x2)Ei(x)  , 

Pi(x)  =  ~  {2e~2z  +  (1  +  2x)e-zEl(x)  +  (1  +  x)  [3£2(x)  -  8P,(2x)] }  , 

.  .  1  (  2P-5P  +  1  ,-6P2+6P  +  l  ,  2  ,  _rr,  ,  /m 

P2(x)=i2  1 - p - e  +( - p2 - +  p?x)e  E^x!p) 


1  9 P3  -  3 P2  —  6P  —  1 

+  (1  +  x)  6(Ei(x)E\(x / P)  -  -Goi(x))  + - y2 - Ex(0x) 


p3(x)  =  (log!  -  l)e  r  -  (1  +  x)  [(logx  -  l)Pi(x)  +  pj2)(x) 

9i (*)  =  \  {2 e~zEl(x)  +  [E2{x)  -  4E1(2x)\ }  , 

,  x  1  fP-1  -Bz  ,-2P2  +  3P+l  P-  l,-XE,//m 
Qli1)  =  2  \  P  e  *  - P2 - p^“x)e  E^xlp) 

II  p3  _  p2  _  op  _  1  ]  'i 

2(E1(x)E1(x/P)  -  — G0i(x))  + - - - £»(/?*)]  j  , 

q3(x)  =  2  {(logx  -  2)e_r  -  [(log i  -  l)Ei(x)  +  Pj2)(x)  }  ■ 

.(78) 

To  calculate  these  functions,  standard  mathematical  software  libraries  have  built-in 
functions  for  the  exponential  integral  and  routines  for  evaluating  the  integrals  (69)  and 
(75)  defining  pj2)  and  Goi-  For  extreme  values  of  the  arguments,  the  appropriate  series 
for  E j  and  G oi  may  be  more  accurate  than  numerical  integration. 


E00x) 
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FIGURE  CAPTIONS 


Figure  1.  Problem  schematic:  Relative  to  the  growing  dendrite  surface,  the  liquid  moves 
steadily  up  and  through  the  interface  as  it  becomes  solid.  Axisymmetry  is  assumed;  gravity 
is  downward.  Far  from  the  dendrite  the  liquid  is  undercooled  and  in  uniform  motion. 

Figure  2.  Parabolic  coordinates:  Lines  of  constant  77  (solid)  and  £  (dashed)  are  shown; 
equal  steps  in  77  and  £  correspond  to  equal  distances  along  the  2-axis. 

Figure  3.  Regions  of  flow:  The  perturbation  solution  applies  within  a  distance  0(G-1) 
from  the  dendrite  tip.  In  the  buoyant  layer  further  downstream  convection  is  nonlinear. 
In  the  far  field  the  flow  is  irrotational  but  not  of  separable  form. 

Figure  4.  Buoyant  flow:  an  example  of  the  instantaneous  flow  field  (velocity  vectors)  near 
the  dendrite  tip,  relative  to  the  solid  material  (the  interface  moves  downward  at  speed  Vr): 
H0  =  1  (S  =  0.6),  P  =  2. 

Figure  5.  Interface  perturbations  A0  (solid)  and  Ai  (dashed)  as  functions  of  H0:  (a) 
P  =  20;  (b)  P  =  2;  (c)  P  =  0.2. 

Figure  6.  Interface  perturbations  Ao  (solid)  and  A\  (dashed)  as  functions  of  P,  for  H0  =  1. 

Figure  7.  Asymptotic  approximations:  For  P  ~2,  Ao  (7a)  and  A\  (7b)  are  shown  (solid), 
with  small-/70  asymptotic  forms  to  0(1)  (as  in  Eq.  56)  (dashed)  and  to  O{H0)  (dotted), 
and  large-f/o  asymptotic  forms  using  one  term  (dashed)  and  two  terms  (see  Eq.  55) 
(dotted). 

Figure  8.  Comparison  with  Ivantsov  solution:  For  one  dimensionless  solution  ( S  =  0.02, 
H0  =  0.004,  P  =  22.8)  the  relation  between  growth  speed  V  (cm/s)  and  tip  radius  R 
(cm)  is  shown,  using  material  constants  for  succinonitrile:  perturbation  theory  (solid), 
and  Ivantsov  solution  (dashed).  The  selection  criterion  a  =  0.0195  is  also  shown  (dotted). 

Figure  9.  Comparison  with  experiments:  Both  perturbation  theory  (solid)  and  the  Ivantsov 
solution  (dashed)  (both  using  the  selection  criterion  a  =  0.0195)  are  shown  with  experi¬ 
mental  data  points  (dots)  from  Huang  and  Glicksman  (1981):  (a)  Peclet  number  Pe  versus 
Stefan  number  5;  (b)  growth  speed  V  (cm/s)  versus  Stefan  number  5;  (c)  tip  radius  R 
(cm)  versus  Stefan  number  5;  (d)  buoyancy  parameter  G  versus  Stefan  number  5. 
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Tip  Solution 


Figure  7b 
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